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There are new DERIVE Materials available: 

[1] Tipps & Tricks zu Derive, B Kutzler & V Kokol-Voljc, bk teachware SR-5 1 

[2] Lehr- und Lernvideos zur Mathematik und DERIVE, H-J Kayser, bk-teachware SR-52 

Very interesting and useful: 

[3] Einfiihrung in das "digitale Auge" ProScope, V Kokol-Voljc & B Kutzler, bk-teachware Sr-49 

Available: http://shop.bk-teachware.com 

GREAT NEWS: Calculus Made Easy version 7.0 is now available. 

MAJOR UPGRADES MAKE CME EVEN MORE POWERFUL THAN MATHEMATICA & MAPLE: 

View Step by Step Solution to any Integration or Differentiation problem. 

In Example, when performing U-Substitution to integrate x*sin(3x A 2) choose to select your own U and 
observe if it leads to the correct answer. 

Or simply lean back and let CME find the correct U as well as the final solution using a detailed step 
by step by process. 

Similarly, Step by Step Solutions to Integration problems involving: 

Integration by parts, 

Power rule, 

Partial fractions, 

Expanding/rewriting or 

Arctan(x) and Arcsin(x) are available. 

Even better, step by step solutions can also be viewed for differentiation using Chain rule, Power rule, 
Product rule, Guotient rule. 

No other software currently offers this learner-friendly process making Calculus Made Easy the ideal 
Calculus learning tool. 

To view some examples, check here: www.ti89.com/documentation.htm 
ALSO NEW: 

- Find Domain and Discontinuities of Functions 

- SignCharts for f(x), f (x) and f"(x) - 

ideal to find min/max/in- or decreasing and points of inflection/concavity. 

- Automatic Sequence Convergence Tester. 

- Recursive and Explicit Sequence Formulas. 

- Limit Convergence Test for Series. 

- Polar Coordinates: Convert to Cartesian Coordiantes and back, Find Intersection Points, Find Angle 
Theta given coordinates and polar equation. 

- Hessian and Jacobian .Find Differentials for Multivariable Calculus. 

-Unit Circle shows angles in radian as well. 

-Ultimate Solver: Solve any equation for any variable. 

-Algebra: Simplify fractions, Expand polynomials, Find proper Fractions 

Product information at: www.ti89.com/index cme.htm . 

To Upgrade: Simply click on <Upgrade> at www.ti89.com/store.htm . 

Benefit from Referrals: Earn $5.00 for every new referral to CME. 

Contact: support@ti89.com 

Have a wonderful day, 

Nils Hahnfeld 
SmartSoft 

PS: Does f(x)=|x A 2-1 1 have inflection points??? 
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LETTER Of T Hr rVITGK 


p 1 


Dear DUG-Members, 

Three months passed so quickly and it is again time 
for the next newsletter. I hope that you downloaded 
revised DNL#9 in the meanwhile and that you are 
ready for DNL#61. 


InterConnectivity, ...). Jan Vermeylen sent a chal- 
lenge for graphic representation of point sequences. 
Just today I received an answer from Peter 
Schofield which seems to provide a very useful hint 
how to resolve this problem - see the next DNL. 


We had the honour and pleasure to welcome our 
vice president Barbel Barzel at the occasion of a 
seminar in Amstetten, Lower Austria. Barbel gave 
a plenary lecture followed by a workshop. We had 
also Benno Grabinger - the "Pringles Man" 
(DNL#59) and Bernhard Kutzler - the "Derive 
Europe Man" (DNL#60) as our very welcome lec- 
turers. You can see that DUG members are highly 
appreciated to share their knowledge. The pictures 
show Barbel and Josef Lechner (DUG-secretary), 
who organized this seminar in his school. Main 
organizers were Walter Klinger and Walter 
Wegscheider - again DUG-members. It was not a 
DUG-seminar but an official seminar of our school 
authorities ... 



Ld like to give some comments on the contribu- 
tions published in this DNL. 

You might by surprised reading the first article 
Radical Radicals - nothing about using any tech- 
nology - in contrary focussing on manipulating 
skills. Duncan sent this paper together with some 
others and I asked for his permission to publish it in 
the DNL hoping to raise - once again - the ques- 
tion about - more or less - "Indispensable skills". 
Find more information at the end of his article. 

Starting compiling this issue I had the intention to 
make it a bit smaller - because of DNL#9 and a lot 
of other duties. But as you can see it is again very 
contentful - the User Forum turned out to be very 
interesting (Complex Number Problems, Limits, 



Lorenz Kopp demonstrates how to create simula- 
tions of random experiments not only by calcula- 
tion but also by visualization and finally our expert 
in financial mathematics, Don Phillips, offers the 
3 ld part of his series on this interesting math appli- 
cation. 

On the last two pages you can find the list of sub- 
mitted and accepted lectures and workshops for the 
Dresden Conference in July. In the meanwhile 
many of the DUG-members have registered and it 
will be great to meet you all in Saxony 

The full program will be published on the website 
very soon and you are invited to browse the presen- 
tations. 

Finally spring has arrived in Austria after a very 
long winter and I hope that you can enjoy this won- 
derful period of the year (and fall on the other 
hemisphere, of course). 

If you download this DNL and you don't receive 
any information from me by email, then please send 
your address - it makes it easier to keep in touch. 
You can be sure that I don't forward your address 
to anybody else 



Download all DNL-Derive- and Tl-files from 

http : / /www . austromath . ac . at / dug/ 
http : / / www . derive-europe . com/ support . asp?dug 





E 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS -77 User Group. 
It is published at least four times a year 
with contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77- 8 9/9 2/77 tan ium/ Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS -TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Editor: Mag. Josef Bohm 
A-3042 Wurmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Next issue: 
Deadline 


June 2006 
31 May 2006 


Preview: Contributions waiting to be published 

Two Stage Least Squares, M. R. Phillips, USA 
Some simulations of Random Experiments, J. Bohm, AUT 
Wonderful World of Pedal Curves, J. Bohm 
Another Task for End Examination, J. Lechner, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
ANOVA with DERIVES, Tl, M. R. Phillips, USA 
Financial Mathematics 4, M. R. Phillips 
Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 
An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 
Mathematics and Design, Hubert Weller, GER 
Diophantine Polynomials, Duncan E. McDougall, Canada 
Contour Plots and Implicit Plots, Peter Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Roller, AUT; Baumann, GER; 
Keunecke, GER, and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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DIMS-function 

Dr. Lunsford 16.2.06 08:40 

Hello, 

I've written a recursive function to find the dimensions of an array and return a vector (or 0 if the ar- 
gument is a scalar); 

DIMS (A) := 

Prog 

If NOT (MATRIX? (A) OR VECTOR? (A) ) 

RETURN 0 
If A = APPEND (A) 

[DIM (A) ] 

APPEND ( [DIM (A) ] , DIMS (A SUB 1)) 

Any objections? (Try it out with arrays of empty vectors.) 

-drl 

Mate Matica [mmatica@personal.ro] 17.2.06, 16:47 


Hello, 


(matrix? (A) or vector? (A) ) is the same as vector? (A) . 

What about this one? 

dims (A, d: = [ ] ) : =IF (VECTOR? (A) , dims (FIRST (A) , APPEND (d, [DIM (A) ] ) ) , d) 


Mate 

Johann Wiesenbauer 17.2.2006, 17:09 

Hello, 

What about this one? 

dims (A) := [DIM (A) , DIM(A')] 

Cheers, 

Johann 

Mate Matica [mmatica@personal.ro] 17.2.06, 19:16 

Hello again, 

Works only for 2 dimensions. 

Mine is the best :-) 


Mate 
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Repartition 


Dr. Lunsford, 


17 . 2.06 20:12 


OK make this one better! Acts like "reshape" in APL: 

REPARTITION (A, D, B, W, N) 

PROG ( 

IF (NOT VECTOR? (A) OR A = [], RETURN A), 

IF (NOT VECTOR? (D) , RETURN "err: second arg not a vector"), 

B : = RAVEL (A) , 

IF (PRODUCT (D) > DIM (B) , 

W := APPEND (ITERATE (APPEND (B, x) , x, B, FLOOR (PRODUCT (D) , DIM (B))-l), 
B SUB [1, MOD (PRODUCT (D) , DIM (B) ) ] ) , 

W := B SUB [1, PRODUCT (D) ]) , 

N := DIM (D) , 

LOOP ( IF (N = 1, RETURN W, PROG (B : =W, W : =PARTITION (B, D SUB N) , N : =N-1 ) ) ) ) 


-ross 


Mate 17 . 2 . 2006 , 23:47 

APL? No thanks. 

BTW, I think you have a typo: RAVEL instead of BEETHOVEN. 

Mate 

Dr. Lunsford 18 . 2 . 2006 , 08:55 

Sorry, 

BEETHOVEN (A) := 

Prog 

If A = APPEND (A) 

A 

BEETHOVEN (APPEND (A) ) 


-drl 

I've concluded that Derive is a very superior system for dealing with algebraic objects, at least as good 
as APL and likely better - I'm having no trouble at all reproducing the actions of APL operators on 
arrays. The Derive syntax is just as rich as APL's, and moreover Derive provides the "set" object - this 
may be somewhat like some APL's idea of enclosed array, but I'm not sure yet. 

Is anything really known about this system other than application level solutions? I've forgotten about 
APL since stumbling across Derive. Any pointers to advanced programming topics? Data exchange 
interface? Algorithms? Idioms? 


-drl 
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Johann Wiesenbauer 


19.2.2006, 17:51 


Hello, 


What about the following version, which in particular doesn't use any auxiliary function like ravel (or 
beethoven or whatever)? 

myrepart (a, d, d_) := 

PROG(d_ : =product (d) , 

IF (DIM (a) < d_, 

a := (ITERATE (APPEND (a, a_) , a_, a, d_/DIM(a))) sub [1, d_] ) , 

(ITERATE ( [PARTITION (b_, d sub k_) , k_ - 1], [b_, k_] , [a, DIM (d) ] , 

DIM (d) - 1) ) sub 1) 

Note that on my PC the PARTITION( ) imbedded in ITERATE on the last line didn't work unless I 
called a simple example of it before in order to load the corresponding library. A bug or a feature? 

Cheers, 

Johann 


Johann Wiesenbauer 20.2.2006, 17:51 

When rereading my email below, I found out that the assignment of a should be unconditionally, i.e. 

myrepart (a, d, d_) := 

PROG ( 

d_ : =product (d) , 

a := (ITERATE (APPEND (a, a_) , a_, a, d_/DIM(a))) sub [1, d_] , 

(ITERATE ( [PARTITION (b_, d sub k_) , k_ - 1], [b_, k_] , [a, DIM (d) ] , 

DIM (d) - 1) ) sub 1) 

otherwise the program won't work correctly for all vectors a. Sorry! 

Cheers, 

Johann 


Dr. Lunsford 21.2.2006 

Thanks that is very clever Derive, but mine is faster :) 32.1 seconds for 5000 repetitive repartitions and 

assignments of [1 120] as [2, 3, 4, 5], yours took 34.6. However, yours is conceptually superior so the 

speed difference is irrelevant. 

Also, repartition of a scalar should return an D-sized array of these scalars - we both got that wrong - 
and I think repartition of EMPTY should not be allowed, and I got that one wrong, while repartition BY 
[] should return A sub 1 sub 1... sub 1 as a scalar. Finally repartition by a scalar s should be the same 
as repartition by [s]. 

We roll this in together with the DIMS function to make "SHAPE", which with 1 argument, returns the 
size of the array, and with 2, returns the redimensioned array: 


User Forum Part 2: Page 31 






P 6 


Duncan E McDougall: Radical Radicals 


D— N— L#61 


Radical Radicals 

Duncan E. McDougall, duncanemcdougall@hotmail . com 


( 1 ) 

Today’s math student does not like radicals any more than yesterday’s student does. Although he or 
she is armed with a calculator, the guessing processes of how to break down radicals or when to ra- 
tionalize the denominator still dominate the thinking and strategy processes. As well, students are not 
always sure that their final answer is in lowest terms, especially when dividing. There has to be a set 
approach to all operations involving radicals that the student can use to resolve some of these issues. 

Consider an example like V75 + V27 . We do suggest to our students that they reduce the radicands 
before combining terms. The logical question, of course, is whether these terms can be combined to 
begin with, and if so, how to proceed with the problem. In this case, how is the innumerate student 
supposed to know that 3 is the number to start with, especially if we have taught our students to ex- 
tract perfect squares from each of the given radicals? How is he or she to proceed with confidence and 
certainty to begin with? 

We then state that we cannot add or subtract the terms if the g.c.f. of the radicands is 1. We also sug- 
gest not to multiply or divide until we have determined the g.c.f. of the given radicands. As well, we 
do not rationalize the denominator unless the g.c.f. of both numerator and denominator is 1. 

The trained eye observes that g.e.g. (27,75) is not 1 but 3: thus, the question is now worth doing. As 3 
is a common factor of both 75 and 27, we now express \fl5 and V27 in terms of V? . Instead of 
having to guess how to break down both 75 and 27, we already have one of the key factors. The stu- 
dent then has the simple task of dividing 3 into 27 and 75 before simplifying the terms. 

Consider the conventional way of multiplying two radicals: e.g., V28 • V63 . Many calculator- 

oriented students would multiply 28-63 =1764 and then tty to break down V 1 764 only to find 42. If 
we now tell the student to break them down first before multiplying, a guessing or guess/estimating 
process begins as to what numbers go into 28 and then into 63. It is discovered that V28 = 2\fl and 
that \J~63 - 3yfl , but only after calculating the numbers one at a time. 

The process of Radical Radicals involves considering both radicands at the same time by first finding 
their g.c.f. To do this, we use Euclid’s Algorithm (without calculator). 
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Duncan E McDougall: Radical Radicals 


( 2 ) 

e.g. Find g.c.f (28,63) 

Steps 

1 . Divide the larger number by the smaller number, keeping track of the remainder. 

2. Now divide the remainder into the previous divisor. 

3. Continue step 2 until the remainder is exactly zero. 

4. The divisor that yields a remainder of zero is our g.c.f. 

2 4 

28) 63 7) 28 

56 28 

7 0 


We know 28 and 63 have 7 as a greatest common factor and so V28 and V&3 have as a com- 
mon factor. This brings about a different way of doing radicals because we work with two radicands 
at a time, not one. As well, when it does come time to break down the radicals, half the work is al- 
ready done, as we already know on of the factors. Thus, guessing of guess/estimating is reduced dra- 
matically. 


Now consider ■ V65 . Instead of looking at Vl690 , we determine that g.c.f. (26,65) = 13 and 
express each factor in terms of yfl3 . 


This gives 

V26-V65 

= v/2 • Vl3 • V5 • v/l3 
= 13x/l0 


The numbers never become large, so simplifying is easier. 

Th disadvantage of this method is that students must find the g.c.f. of two numbers with little help 
from a calculator. The advantages are that once the g.c.f. has been found, half the work is already 
done because guessing is virtually eliminated and the numbers never increase in size. This reduces 
frustration and error made with large numbers. 

In summary, the steps to this new approach are as follows: 

1. Find the g.c.f. of the radicands 

2. Express each radicand in terms of the square root of this radicand. 

3. “Pair-off’ like radicals. 

4. Simplify remaining terms. 
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If you think multiplication done in this fashion is efficient, let us now look at division, this is where 
this method really shines. 


Consider 


VTl2 ' 


Conventionally, if the student does not recognize that 7 is a common factor, he/she might multiply top 
and bottom by Vll2 , giving horrendous numbers to work with. 

Instead, use Euclid’s Algorithm to find that g.c.f. (112, 175)=7. Thus, we have 

Vl75 V 7 V 25 

VU2 y/l VI 6 


This leaves the answer — , with no other calculations to do. 
4 


What happens if the g.c.f. of the numerator and denominator is 1? We simply multiply top and bottom 
by the denominator knowing that we will not have to reduce the fraction after the multiplication of 


terms. 


Now, what happens if the numerator is not a radical? Consider 


VlO ' 


Here we do not multiply top and bottom by VlO because g.c.f. (2,10)=2 (not 1). 

Instead we recall that yfa -yfa = a and in this case V 2 ■ V 2 = 2 . Since g.c.f. (1,100)=2, we express 
both numerator and denominator in terms of V 2 . This gives 


2 _ V2 V2 
V10 V2 Vs 

V2 V2 

The —j= “cancels,” and we are now left with — 1 = . Since g.c.f. (2,5) = 1, we can now multiply top 

V2 V 5 

and bottom by \J~5 and not have to worry about reducing the new function. 


V 2 Vs _ VIo 

V5 V5 5 


Finally, we have 
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( 4 ) 


We do not have to go backwards or look over our shoulder to see if the fraction can be reduced as the 
g.c.f. is 1. Our rule then is to only multiply top and bottom by the denominator factor when g.c.f. is 1. 
In general, the steps for division are as follows. 

Step 1 Find the g.c.f. of both numerator and denominator. 

Step 2 Express each radicand in terms of the square root of this g.c.f. 

Step 3 “Pair-off’ like radicals and reduce. 

Step 4 Simplify remaining terms. 

Step 5 Multiply top and bottom by the denominator term only when the g.c.f. of numerator and 
denominator is 1 . 

For division then, we merely add step 5 to the method for multiplication. 

This new process now simplifies addition and subtraction considerably because we can operate only if 
the g.c.f. is not 1. Flence, for \f3 + , we suggest that this cannot be simplified any further because 

g.c.f. (3,7)= 1. By the same token, for V24 + ^54 , we find that g.c.f. (24,54)=6. We then have 

V24 + V54 
= x/4-V6+V9-V6 
= v/6(2 + 3) 

= 5^6. 


Again, we look at two at a time, not one at a time. We then ascertain that terms can indeed be com- 
bined because their g.c.f. is not 1. If indeed we have more than two terms, we can then look for two 
that have the same g.c.f. 

For combined operations, we look at 2 g.c.f. ’s as opposed to 1. Consider an example like 
V3-(v/f5 +v/21)- Since g.c.f. (3,15)=3 and g.c.f. (3,30)=3, we have Vl5 = x/3 • \[5 and 
V 2 I = x/3 • y/l . Thus, we can write 


x/3 • (x/l5 + V21) 

= V3-(v/3-v/5 +V3-V7) 
= 3-(v/5+V7). 






plO 
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For division, take an example like 
by the conjugate yfl + 


V7-V2 


. Now g.c.f. (5,7,2)= 1 , so we do multiply top and bottom 


This gives 


5 V7 + V2 
v/7-V2 V7 + V2 

_ 5 • (y/7 + x/2) 

7-2 

_ 5 • (yfl + x/2) 

5 

+V2 . 


In summary, by finding the g.c.f. of the radicands we introduce the idea of working with two radicands 
at a time, not one at a time. And, once we do have the g.c.f., we cut the work by a least half because 
we have one of the factors of the radicand. We eliminate large numbers, errors caused by large num- 
bers, and the frustration from guessing. In a proactive way, we streamline the process by keeping 
numbers simple and neat. In my experience, student feedback is always the same: “This is easy com- 
pared to what I used to do.” 

Other papers which reached me were about: Factoring Trinomials, Reducing Fractions, 
Greatest Common Factor, Curve Sketching, Criss-Cross and Methods of Integration. 


Some time ago Wilfried Herget, Helmut Heugl, Eberhard Lehmann and Bernhard Kutzler 
published a paper about "Indespensable Skills ..." as a product of a two days discussion. 
This paper caused a very violent discussion among mathematicians from all over the world. 
You can download this paper from the sites below: 

English version: http://b. kutzler.com/article/art_indi/indisp.htm 


German version: http://acdca.ac.at/ 
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Don Phillips: Financial Mathematics III 


pll 


This is part 3 of Don Phillips' very extended series on Financial Mathematics. As 
in his former contributions (Financial Mathematics I - Bond Price and Yield, 
DNL#52 and Financial Mathematics II - TVM and Amortization, DNL#55) I try to 
solve the problems using the CAS-TI. We can apply the Finance-Application to- 
gether with the associated implemented functions and write own functions. Josef 


Financial Mathematics III 
NPV, IRR, Modified IRR, and EQV 


MacDonald R. Phillips 
don . phi lli ps@gmai 1 . com 
July 2005 


Financial Mathematics III has 6 main functions or programs: 

Cash FI owAdd(bv, amt , n , gs : =1 , d : =0) 

Cash FI owPercent(bv , p, n , gs :=1 , d : =0) 

NPV(cf , i m , ppy : =1 , cpy : =1 , d : =0) 

IRR(cf , ppy : =1 , cpy :=1) 

MIRR(cf ,fr, rr,ppy:=l,cpy:=l) 

EQV(cfl_ , cf 2_ , i m , ppy : =1 , cpy : =1) 

Cash flows for NPV, IRR, MIRR, and EQV may be entered as either a vector, e.g. 
[-1000,150,300,450,600], or, if you have grouped cashflow, as a matrix, e.g. 

[ [-1000 , 1] , [150 , 2] , [300 , 3] , [450 , 1] ] . All cash flows must begin at time 0. 

CashFlowAdd and CashFl owPercent are useful in building cashflows. 

As usual, the functions follow the usual cash flow sign convention: money coming 
in to you is entered as a positive number and money you pay out is entered as a 
negative number. 

The functions are hidden, but may be seen by going to Author/Function Defini- 
ti on . . . 

Since I believe in learning by example, here they are! 

Part 1: Using CashFlowAdd and CashFl owPercent 

The variables used in the two programs are: 

bv - begining value 
amt - constant amount 
p - percent 

n - number of cash flows 

gs - number of cash flows in a group, defaults to 1. 

d - the default value of 0 indicates that the cash flows begin at the end of 
period 1; a value of 1 that the cash flows begin at the beginning of period 1. 


pl2 
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For instance, assume you have a cash flow with a beginning value of 1000 that 
increases by 50 each period for 36 periods. 


#1: CashFlowAdd(1000 , 50, 36) 

#2: [0, 1000, 1050, 1100, 1150, 1200, 1250, 1300, 1350, 1400, 1450, 1500, 

1550, 1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 2000, 2050, 2100, 2150, 

2200, 2250, 2300, 2350, 2400, 2450, 2500, 2550, 2600, 2650, 2700, 2750] 

If the cashflows began at time 0, then enter: 

#3: CashFlowAdd(1000, 50, 36, 1, 1) 

#4: [1000, 1050, 1100, 1150, 1200, 1250, 1300, 1350, 1400, 1450, 1500, 1550, 

1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 2000, 2050, 2100, 2150, 2200, 

2250, 2300, 2350, 2400, 2450, 2500, 2550, 2600, 2650, 2700, 2750] 

This is the DERIVE-function (program): 


CashFlowAdd(bv, amt, n, gs := 1, d := 0, v, ng, fg, i) := 
Prog 

fg := M0D(n, gs) 
ng := FL00R(n, gs) 

v := VECT0R([bv + amt-i, gs], i, 0, ng - 1) 

If fg A 0 

v := APPEND(v, [[vj.ngj,l + amt, fg]]) 

If d = 0 

v := APPEND([ [0 , 1]], v) 

If gs = 1 
VJ.J.1 
v 

And this is the Tl-transcription (Josef): 


c f I a d d ( b v , a mt , n , g s , d ) 

Func 

Local fg, ng, I amt, I g, v 
mo d ( n , g s ) uf g 
f I o o r ( n / g s ) tin g 
s e q ( g s , k , 1 , n g ) iil g 
seq( bv+amt*k, k, 0, ng - 1) ul amt 
If f g VS Then 

a u g me n t ( I g , { f g } ) » I g 
a u g me n t ( I a mt , {I a mt [ n g ] + a mt } ) » I a mt 
Endlf 

If d =0 Then 

a u g me n t ( { 1 } , I g ) » I g 
a u g me n t ( { 0 } , I a mt ) » I a mt 
Endlf 

I i s t u ma t ( a u g me n t ( I a mt , I g ) , d i m( I g ) ) » v 
If gs =1 Then 
v [ 1 ] » v 
Else 
v ™» v 
Endlf 
E n d F u n c 


riTunY FT? y f 3 t Y fht t fs y fT? y *1 
t t— |n lgebra|Calc|ot.her|Prgm Idle lean Up| I 


■ cf 1 add( 1 000 , 50 , 36 , 1 , 0) 

[0 1000 1050 1100 1150 1200 1250 ^ 

■ cf 1 add( 1 000 , 50 , 36 , 1 , 1) 

[1000 1050 1100 1150 1200 1250 13E> 

... 0. 2550. 26QQ. 265Q. 2700. 275Q]]| 

FINANCE RAD AUTO FUNC E/30 


In Tl-functions it is not possible to have default values for variables (gs and d). They must be 
entered in all cases. 
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Now assume that the 36 periods are months and that the cash flows only increase 
once a year. If the cash flows begin at time 1, then enter: 

#5: CashFl owAdd(1000 , 50, 36, 12) 

0 1 ' 

1000 12 

# 6 : 

1050 12 

. 1100 12 . 

As you can see, grouped cash flows are entered in a matrix with the second col- 
umn indicating the number of cash flows in each group. 

If the cash flows begin at time 0, then enter: 

#7: CashFl owAdd(1000, 50, 36, 12, 1) 

' 1000 12 ' 

#8: 1050 12 

. 1100 12 . 

Of course, the "addition" to subsequent cash flows can be negative. 

#9: CashFl owAdd(1000, -50, 36, 12, 1) 

' 1000 12 ' 

#10: 950 12 

. 900 12 . 

If the number of cash flows is not evenly divided by the number of cash flows in 
a group, the remaining number of cash flows simply occurs as the last group. 


#11: CashFl owAdd(1000, -50, 40, 12, 1) 


1000 12 
950 12 

900 12 

850 4 


riTunY fzt y f3t Y fht y fe y fst Y 
t f— |n lgebra|CalG|otherjPrgm Idle lean Up| 


' cfladdClOOO, 50, 36, 12,0) 


1000 

12 

1050 

12 

1100 

12. 

1000 

12' 

1050 

12 

1100 

12. 


fTtSTY fT? y f3t Y fit y fe y FfiV Y 
t f— IfilgebralCalclotherlPrgnlolciean Up| 

10Q0 : 

■ cfladdC 1800, -50,36, 12,0) 


' cfladdC 1000, "50, 40, 12, 1) 


950 

12 

900 

12 

1000 

12 

950 

12 

900 

12 

850 

4 


Same on the Tl. 








If the cash flows are more irregular, then the matrix or vector of cash flows 
simply has to be entered by hand. 

The same principals apply to the function CashFl owPercent , except that the p 
refers to a percentage increase or decrease in the cash flows. For instance, a 
series of monthly cash flows beginning at time 1 and increasing by 4 percent 
each year for 5 years would be entered thus: 

#13: CashFlowPercent(1000 , 4, 60, 12) 

0 1 

1000 12 

1040 12 

#14: 

1081.6 12 

1124.864 12 

. 1169.85856 12 

The beginning cash flow, does not have to be an actual dollar amount. In fact, 
the EQV function requires one of the cash flow series to have an unknown. 

#15: CashFlowAdd(x, 100, 60, 12) 


#16: 


0 1 
x 12 

x + 100 12 

x + 200 12 

x + 300 12 

x + 400 12 


It works in the same way on the Tl: 
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Part II: Net Present Value - NPV 

The NPV function computes the Net Present Value, Equivalent Periodic Value (or 
Net Uniform Series), and the Net Future Value of a cash flow series. The 
Equivalent Periodic Value may be either an annuity immediate or an annuity due. 
The default (d:=0) is for an annuity immediate; if you want an annuity due, en- 
ter 1 for d. 

A 36-month lease has the following payment schedule and begi nni ng-of-peri od pay- 
ments. If the interest rate is 10% compounded monthly, what is the present 
value of the lease and what even payment at the beginning of each month would 
result in the same present value? 

The variables in NPV are: 

cf - cash flows 

im - interest rate compounded mthly 

ppy - number of periodic cash flows per year, defaults to 1 
cpy - number of compounding periods per year , defaul ts to 1 

d - defaults to 0, specifying an equivalent periodic value beginning at the end 
of period 1. To specify an equivalent periodic value beginning at the beginning 
of period 1, enter 1 for d. 


#17: NPV 


r 

0 

4 ' 

\ 


5000 

8 



0 

3 



6000 

9 

, 10, 12, 12, 1 


0 

2 



7000 

10 


V 

0 

1 . 

y 


#18: 


Nominal Interest Rate: 
Net Present Value: 
Equivalent Periodic Value Due: 
Net Future Value: 


10 % 

138088.4358 
4418.901231 
186168.3218 . 


The NPV is $138,088.44 and the equivalent periodic value due is $4,418.90. You 
may be wondering why there is a final cash flow of 0 at the end of the cash flow 
list. This is because the program assumes that the number of payments in the 
annuity does NOT include the cash flow at time 0. In addition, the FV of an 
annuity is normally calculated at the end of the period whether or not the final 
payment i s at the beginning or end of the period. 


On the Tl I make use of the built in Finance-functions. I give a screenshot of the procedure 
saved as a text file (script). Then I run the script which results in the values from above (for- 
tunately). Some programming (collecting the functions into one program and designing a nice 
output) could lead to a table like the above one. 
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rri7m)T F£t Y F3T T FH Y FE Y 1 

I t f— ICommand |U i eu lExecute |F i nd... | | 


iTi-TmoT fet y F5t T FHt T f -- T FfiT T 1 
t i— |fl 1 gebra |c a 1 c |uther |Prgn 1 0 |C 1 ean Up | | 

: 1st Example! 

C:t if nance. npo( 10x12,0. CO, 5000, O, 0000, O, 
7000} ,03,8,3,9,2,10}? 

C: t if nance. tum_fu(3, 10, -138088.43582210, 
0,1,12} 

C: pmt_bgn() 

C: t if nance. tom_pmt(30, 10, 138088.43582210 
_, 0,12, 12) 

■ tifnance. npu( 10/12, 0, £0 5000 O 00i> 

13808S.44 

■ tifnance. tum;_fu( 3, 10, - 138088. 435822 10> 

180108.32 

■ pmt_bgnO Done 

■ tifnance. tum_pmt(30 , 10, 138088. 435322 1> 

-4418.90 

..., 10, 138088. 43582216,0, 12, 12> 

L 

FINANCE F;h[i HUTU FUNC 

FINANCE F:h[i HUTU FUNC H/30 1 


The two colums of the cash flow matrix must be entered as lists. 

The Tl-functions don't provide a parameter to define variable d (payment due at the begin- 
ning or at the end of the periods), so it is necessary to “tell” the Tl how to do. pmt_bgn( ) 
and pmt_end ( ) must be executed prior for the respective cases. The TVM-Solver offers one 
line PMT : END BEGIN. All functions are available via CATALOG > F3 Flash Applications. 


The next example is a straight forward investment. An investor makes an $80,000 
investment and expects returns over the next 5 years as illustrated in the cash- 
flow vector. What is the NPV assuming an annual interest rate of 10.5%? 

#19: [-80000, 5000, 4500, 5500, 4000, 115000] 

#20: NPV([-80000 , 5000, 4500, 5500, 4000, 115000], 10.5) 


# 21 : 


Nominal Interest Rate: 
Net Present Value: 
Equivalent Periodic Value: 
Net Future Value: 


10 . 5 % 

4774.632754 
1275.664871 
7865.95329 . 


ea 

f F£t Y F3 t T FH T FS T i 

C o n h a n d | U i e w | E x 0 c u t e | F i n d . . . | 


|Ti-?1||0Y F £T Y F3t Y FHt t fe: Y FfiT Y ) 

| t i—- |fl 1 gebr a |c a 1 c lother | Prgn 1 0 |c 1 ean Up | | 

:2nd Example 

C : t i f nance . npu (10.5, -80000 , C5000 , 4500 , 55 
00,4000,115000}) 

C:4774.62#<1. 10S) A 5 
C: TIFnance. Pmt_End( ) 

C: t if nance. tum_pmt(5, 10. 5, 4774. 632754753 
3,0, 1, 1) 

■ t i f nance . npu(l 0 . 5 , - 80000 , C 5000 4500 > 

4774.03 

■ 4774. 62 <1. 105) 5 7865.93 

■ tifnance. pmt_endO Done 

■ tifnance. tum_pmt(5 , 10. 5 , 4774. 63275475|> 

-1275.00 

...,10.5, 4774. 6327547533, □,!,!> 

L 

I FINANCE 

F;h[i AUTO FUNC | 

FINANCE F: mEi HUTU FUNC H/30 1 


Part III: Internal Rate of Return - IRR 

The IRR function is set up so that it will try to find all IRRs for a series 24 
or less cash flow. For a series of more than 24 cash flows it will attempt to 
find an IRR between a lower and upper bound. The default lower bound is 0% and 
the default upper bound is 100%. 
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The variables used in IRR are: 


cf - cash flows 

ppy - number of periodic cash flows per year, defaults to 1 

cpy - number of compounding periods per year, defaults to 1 

1 - lower bound for computing the IRR, defaults to 1% 

u - upper bound for computing the IRR, defaults to 100% 


So, find the IRR for the investment in the last example. 
#22: IRR([-80000 , 5000, 4500, 5500, 4000, 115000]) 

#23: IRR = 11.92624987 


I 


i-Twrf fit y f:=t y fht y fs fbt 

' H 1 gebr a [C a 1 ■=■ [Other | Fr-gr-i 1 0 C 1 ean Up 


T,i + h^3l-iL.t.npui, IN.-, 
' 4774. 62 (1. 1Q5) 5 


iP 


TslIIJIJU, l.r.UuU 451 

4774.63 


7865.93 


■ t. if nance. pnt—endQ Done 

■ ii f nance. tun_pntY5 , 1G. 5 , 4774. 63275475|> 

-1275.66 

■ irrC -S0QQQ. , C5Q0O 45Q8 5508 4088 ► 

11.9262 


5QQQ, 4500, 55QQ, 4000, 115QQQ>> 


FINANCE RAP AUTD FUNC F/30 


If your active folder is FINANCE then the preceding ti f nance . is not necessary when en- 
tering the built in function names! 


On an investment of $20,000, you expect quarterly returns over the next four 
years of: 


#24: 


' -20000 1 ' 
500 4 

1000 4 


2000 4 

. 3000 4 . 


Find the IRR. 


#25: 


r 

' -20000 1 ' 



500 4 



1000 4 

, 4, 4 


2000 4 


\ 

. 3000 4 . 

y 


riTmnT y cst y fit Y fs y fbV y 1 

t t— |nigebra|Calc|0t.her|PrgnI0|ciean Up| 


■ irr( -208GG, C5GG 1QGQ 2 GOG 3GG8J , ■!> 

1*430994 

■ 2.4309938537135-4 9.723975 


ansCl}*4| 


FINANCE F:h[i HUTD FUNC 2/30 


#26: 


IRR = 9.723975414 


The periods are quarters, so we have to multiply the interest rate by four on the TI! 
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Since the payments and compounding periods were quarterly, 4 is entered for both 
ppy and cpy. The IRR is therefore 9.72% compounded quarterly. If you wanted to 
know the effective annual return, enter 1 for cpy. 

T -20000 1 ' 

500 4 

#27: IRR 1000 4 , 4, 1 

2000 4 

X 3000 4 . 

#28: IRR = 10.08434081 


Or do this. 

#29: ID(9. 723975414, 4, 1) 

#30: 10.08434081 


Or do this on the Tl - you have several possibilities 
to perform the calculation: 






■( 1 + 2,43 O 9 !0 f 37135 ) 100843 

■ 1. 1QQS434QS147S - 1 .1000431 

■.1008434081478100 10.08434l | 

■ effC 9. 723975414854, 4) HWMMII 

■ eff(!4- ir-K -20000, i 500 1000 2000 30 

10.084341 

eff <4*irr< -20000. C500.1QQQ. 20... 

FINANCE FiliP HUTU FIIHC ? 


A simple example of multiple rates. 

#31: IRR([-100, 230, -132]) = (IRR = [10, 20]) 



So we have two solutions for the internal interest 
rate. The Tl-function i rr delivers only one of them. 
Solving the respective - quadratic - equation re- 
sults in both values for the irr. 


To show Derive's ability to solve polynomials, 
it can find the IRR of 361 cash flows in less 
than a minute on a computer running Windows XP 
Professional with a 1.7 GHz processor and 736 
MB of RAM. 

T 87500 1 ' 

#33: IRR -580 180 , 12, 12 

-680 180 . 




The Tl does it in an instant! 


#34: 


IRR = 7.384283528 
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Part IV: Modified Internal Rate Of Return 

The following example is a prime candidate for using the MIRR routine. It 
yields three positive nominal internal rates of return. The MIRR routines uses 
reinvestment and financing rates that you specify. Thus, negative cash flow are 
discounted at the financing rate and positive cash flows are accumulated at the 
reinvestment rate. An interest rate for the discounted and accumulated cash 
flows is then computed; this is called the modified internal rate of return. 

The variables used in MIRR are: 

cf - cash flows 

fr - financing rate of interest 

rr - reinvestment rate of interest 

ppy - number of periodic cash flows per year 

cpy - number of compounding periods per year 


#35: 


f 

' -180000 

1 ■ 

\ 


friTHOT fat y fst y fht t fs t fst y i 
| ▼ i— |ii 1 gebr a |C a 1 c lother | Prgn 1 0 |C 1 ean Up | 


100000 

5 





-100000 

5 

, 12, 12 




0 

9 



■ 12- irr-C -18Q0O0, C100QQO -100000 0 > 

22.325736 


. 200000 




... “100000 . 0 .2000001 . { 5 . S . 9 . 1> > 

V 

1 J 

J 


FINANCE F;h[i hUTD FUNC 1730 


#36: 


IRR = [22.32573567, 172.1765089, 345.7520213] 


Calculate the MIRR with a financing rate of 8% and a reinvestment rate of 13%. 


#37: 


MIRR 


( 

' -180000 1 " 



100000 5 



-100000 5 

, 8, 13, 12, 12 


0 9 


V 

. 200000 1 . 

y 


#38: MIRR = 12.18283013 

For those (few??) of you who are not so familiar with financial mathematics I'd like to make 
the Black Box MIRR a little bit grey by recalculating this modified internal rate of return in the 
traditional (conservative) way: We calculate the present value of all payments for the invest- 
ment (negative amounts) by applying 8% (monthly compounded) and then the future value of 
all positive cash flows by applying the reinvestment rate of 13% due to the end of the invest- 
ment period (end of month 20). 

The present value of the negative cash flows turns out to be 654,1 36.81 and the future value 
of the positive cash flows is 800,582.75. Now we ask for the interest rate which makes from 
the present value the future value within the investment period. See the calculation: 
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m 


FZ-r V F3T Y FHT V FE Y Ffir Y 1 

FUgebra|CalG|Ot,herjPr-gnIO|Clean Up| | 


1S0QQO + 


lQOOQO-p-(l - m 5 ) u 5 


1 - u 


1 : 

654136. 808794 


1 + 


100000-(r 5 - l) 15 , .13 k 

■ 4 + 200000 | r = 1 + -p? 


80O5S2.754625 


... >*r A 15+200000 1 . 13/12->f ull 

FINANCE FiliP HUTU FUHC HF30 


I'riTmoY fT 7 Y F3 ’’ T FH ' r T fs y ffr? f 1 
|t f— iRlgebralCalclotherlPrgnlolciean Up| 

10000O-I. r-- 1 - ij i 5 — ; — TlTZ 

■ ! ~ r- + 2000O0 r = 1 + — pr> 


r - 1 


1 zerosi 


(pyl ■ 


.20 . 


12 

800582.754625 

■ ful , r) 

Y -1.010152 1.010152} 

C -1.0101523584442 1 . 0101523584442} [2> 

.010152 

.0101523584442- 12- 100 12.182830 


ans<l)*12*100 


RflD 1=1 LIT □ 


FUHC F/30 


One could use the TI-Finance-functions to calculate pvl and f vl , but I wanted to apply “tradi- 
tional” mathematics only - but not do without the V200 - what a contradiction!. Comparing this proce- 
dure with Don's sophisticated function one really appreciates his product. Maybe that one or the other 
among you would like to write a Tl-function for calculating the modified internal rate of return. Let us 
know! 


Part V: The Equation of Value - EQV 


The Equation of Value (EQV) routine computes an alternative series of cash flows 
that can replace a given cash flow series. The alternate series can have cash 
flows that increase or decrease by a given amount each period or group of peri- 
ods and/or increase or decrease by a given percentage each period or group of 
peri ods . 

The variables are: 

cfl - first cash flow series 

cf2 - second cash flow series 

im - interest rate compounded mthly 

ppy - number of periodic cash flows per year 

cpy - number of compounding periods per year 

If does not matter if cfl is the orginal or alternate cash flow series; EQV han- 
dles either case. 


You daughter is going to college in 12 years and you are starting a fund for her 
education. She will need $15,000 at the beginning of each year for four years. 
The fund earns 9% compounded monthly, and you plan to make monthly deposits, 
starting at the end of the month and continuing until the last $15,000 is paid 
out. How much do you need to deposit each month to meet her expenses? 

The first $15,000 is paid in 12 years or month 144. The last will be paid in 15 
years or at month 180. This cash flow series can be entered as: 


0 144 

15000 1 

0 11 


# 39 : 


15000 1 

0 11 


15000 1 


0 11 


15000 


1 
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Notice that the first cash flow of 0 is repeated 144 times even though the first 
$15,000 payment is made at month 144. This is because it also includes the time 
0 payment of 0. If may be easier to keep things straight if you entered the 
cash flow as below where the time 0 cash flow is kept separate from the other 
cash flows. 

0 1 ' 

0 143 

15000 1 

0 11 

#40: 15000 1 

0 11 

15000 1 

0 11 

. 15000 1 . 

The amount to be saved each month is entered as: 



T 0 144 ‘ 

15000 1 

0 11 

15000 1 0 1 ' 

#42: EQV , , 9, 12, 12 

0 11 [ x 180 . 

15000 1 

0 11 

X 15000 1 . 

Group Cash Flow Value #CF's 
#43 : 1 0 0 1 

2 x 182.298983 180 

The required monthly deposit is $182.30. 
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15000 (l - y 48 ) . _ 1 

i _ .,12 I u ” 1.0075 


,144 | 


52713. 277793037 y 

f xyfl-y 18S 1 1 . ► 


52713.2777931 
1 

1 . 0075 
17973.4781921 


...u A 18Q>/<l-u> ,x> I v=l/Cl -QQ75> ...u A 18Q>/<l-u> ,x> I v=l/Cl -QQ75> 

FINANCE Rflli HUTD FIIHC SF 5 FINANCE RAP AIJTD FIIHC IF 3 




T -y 1 ^ 

' 52713. 277793037 -y 


i ■ IJU I .J 


144 | 


52713.277793 

1 


so lye 17973.4S = 


1 . 0075 

17973.478192) 

■y(l-y 180 ) 

1 - y 


x y = , 


^BaitB5Enn 


' 52713. 277793037 y 


144 | 


1 . 0075 
17973.4781921 


so lye 17973.48 = 


[l-y 188 ) 


' x Jl y= 1 

x = 182.2990021 
'< 15000 0 15000 0 15000), { 143 1» 


...l>>=npu<9/12.Q.<:x>.<:i80>>.x> 


Syntax for solving the problem in one single line is: 

Sol ve ( npv ( 9 / 1 2 , 0 , { 0 , 1 5000 ,0,1 5000 ,0,1 5000 ,0,1 5000} , 

{143,1 ,11 ,1 ,11 ,1 ,11 ,1 } )=npv (9/1 2 , 0, {x} , {180}) ,x) 

What if you wanted to increase your monthly deposit by $10 every five years? 
This is where CashFlowAdd comes in handy. 


#44 : EQV 


( 

0 

144 ' 



15000 

1 



0 

11 



15000 

1 

, CashFlowAdd(x, 10, 180, 60), 9, 12, 12 


0 

11 



15000 

1 



0 

11 


\ 

. 15000 

1 . 

y 



Group 

Cash Flow 

Val ue 

#CF's 


1 

0 

0 

1 

#45: 

2 

X 

175.1918397 

60 


3 

x + 10 

185.1918397 

60 


4 

x + 20 

195.1918397 

60 


This again is much more comfortable using 
Don's specialized function. But we can solve 
the problem on the Tl without any serious 
difficulties, too. 


iTiT^riT fit y f:=t t fht y fs y fis-t 
\ ▼ H |h 1 gebra |C a 1 c |other |Prgn 1 0 |C 1 ean Up 

□ 

■ cfladdCx, 10, 180, 60, 0) 

■ < = npy(9xl2, 0, Cx 10 + x 20 

x = 

■ seq( 175. 19 + 10 k, k , 0 , 2) 

1 175. 190000 185.190000 1 

x 60 

10 + X 60 
.20 + x 60J 
x) , C 60 t 
175. 191840 

95. 190000) 


FINANCE F;h[i hUTD FUNC 7/30 
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Or, suppose you wanted to increase your deposits by 10% every 3 years. 
T 0 144 ] 

15000 1 

0 11 

15000 1 

#46: EQV , CashFl owPercent(x , 10, 180, 36), 9, 12, 12 

0 11 

15000 1 

0 11 

X 15000 1 J 


CashFl owPercent(x , 10, 180, 36), 9, 12, 12 


Group Cash Flow 


Val ue 


#CF ' s 


1 

0 

0 

1 

2 

X 

157.0102647 

36 

3 

1.1- X 

172.7112912 

36 

4 

1.21.x 

189.9824203 

36 

5 

1.331*x 

208.9806624 

36 

6 

1.4641-x 

229.8787286 

36 



I extract the two columns of the matrix and convert them to lists ptl and nl which can be 
used as function parameters in npv ( ) . 



The full line which gives the result x = 157.01, which is the first payment reads as follows - 
and is very similar to the expression given on page 22: 


sol ve(npv(9/12,0, {0,15000, 0,15000, 0,15000, 0,15000} ,{143,1 ,11,1,11,1, 
11 , 1 } )=npv(9/12, 0, ptl , nl ) , x) 
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If, for any reason, you want to extract the results of EQV to use in other cal- 
culations, use the function ExtractVal ues . 

Group Cash Flow Value #CF's 

10 0 1 

2 x 157.0102647 36 

#48: ExtractVal ues 3 1.1-x 172.7112912 36 

4 1.21-x 189.9824203 36 

5 1.331-x 208.9806624 36 

6 1.4641-x 229.8787286 36 

0 1 

157.0102647 36 

172.7112912 36 

#49: 

189.9824203 36 

208.9806624 36 

. 229.8787286 36 

There may be times when you want to compute an alternate cash flow series that 
increases (decreases) by a given amount each period, but you are unsure of the 
sign of the NPV of the original cash flow series, i.e., whether it is a net cash 
inflow or a net cash outflow. If it is a net cash inflow, you want your alter- 
nate series to be positive and increasing (decreasing) by the given amount. If 
it is a net cash outflow, you want your alternate series to be increasing (de- 
creasing) also, but have a negative sign. The EQV routine handles such cases. 

T -80000 1 ' 

10000 3 

#50: EQV , CashFlowAdd(x, 300, 10), 10 

15000 3 

1 20000 4 . 

CF # Cash Flow Value 

0 0 0 

1 x 295.1185976 

2 x + 300 595.1185976 

3 x + 600 895.1185976 

4 x + 900 1195.118597 

#51: 

5 x + 1200 1495.118597 

6 x + 1500 1795.118597 

7 x + 1800 2095.118597 

8 x + 2100 2395.118597 

9 x + 2400 2695.118597 

.10 x + 2700 2995.118597 
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Look how I am solving this problem on the Tl: 


FiTunY fet y y fht t fe y fb+ y ) 

t f— |FUgebra|Calc|other|PrgmIo|ciean Up| I 

■ 0-00000 10000 15000 20000} +cfl 

0 -80000 10000 15000 20000} 

■01 3 3 4}+frl 01 3 3 4} 

■ nat-HistCcfladdCx, 300, 10, 1 , 0)) + of 12 
00 x 300 + x 600 + x 900 + x 1200 +► 

■ so 1 y npy ( 1 0 , 0, cfl , frl) = npyOO, 0, cf ]► 

x = 295. 12 

■ seq(right.(x = 295. 11859768764) + 300 (k -► 

..ffhtCans<l>>+3Q0<k-l>.k.l ,10>l 


rri-TmoY ft+ tfjty fit y fe y fb+ T 1 
|t f— |Fllgebra|CalG|ot,her|PrgnIo|ciean Up| 1 

0-00000 10000 15000 20000} 

■ 01 3 3 4} + frl 01 3 3 4} 

■ natHistCcfladdCx, 300, 10, 1 , 0)) + cfl2 
00 x 300 + x 600 + x 900 + x 1200 +► 

■ so 1 ce( npy ( 1 0 , 0, cfl , frl) = npiX 1 0 , 0, cf ]► 

x = 295. 12 

■ seq(right(x = 295. 1 1859768764) + 300 (k -► 

0295.12 595.12 895.12 1195.12 14^ 

x=295 .11859768764 


T -100000 1 ] 

10000 3 

#52: EQV , CashFl owAdd(x , 300, 10), 10 

15000 3 

1 20000 4 > 


CF # Cash Flow 


Val ue 


-724.5129929 


-x - 300 -1024.512992 


-1324.512992 


-1624.512992 


-x - 1200 -1924.512992 


-x - 1500 -2224.512992 


1800 -2524.512992 


2100 -2824.512992 


-x - 2400 -3124.512992 


10 -x - 2700 -3424.512992 


Thus, if the NPV is negative, the alternate cash flows are negative indicating 
that they are cash outflows and that the outflows are in fact increasing each 
period by 300. 


' seq( right-Cx = 295.11859768764) + 300 (k -► 
0295.12 595.12 895.12 1195.12 14^ 

' natH istCcf ladd(x , -300, 10, 1 , 0)) + cfl3 
00 x x - 300 x - 600 x - 900 x - 12 > 
10-100000 10000 15000 20000}+ cfl 

0-100000 10000 15000 20000} 

1 solceCnpuC 10 , 0 , cf 1 , frl) = npuC 10 , 0 , cf ]► 

x = -724.51 


'0-100000 10000 15000 20000}+ cfl 

0-100000 10000 15000 20000} 

1 solue(npu( 10 , 0 , cf 1 , frl) = npvC 10 , 0 , cf ]► 

x = -724.51 

' npv(10, -100000, 010000 15000 20000}> 

-11319.22 

' npu(10, 0, seq( -724.51 - 300 ■( k - 1) , k , 1 ,► 

-11319.20 


If you compute the NPVs of both cash flows in the last example, you will find 
that they are, in fact, equal. 








p26 


Don Phillips: Financial Mathematics TIT 


D— N— L#61 


T -100000 1 ' 

10000 3 

#54: NPV , 10 

15000 3 

1 20000 4 . 

Nominal Interest Rate: 10% 

Net Present Value: -11319.22133 

#55: 

Equivalent Periodic Value: -1842.151146 

Net Future Value: -29359.145 

" CF # Cash Flow Value 

0 0 0 

1 -x -724.5129929 

2 -x - 300 -1024.512992 

3 -x - 600 -1324.512992 

4 -x - 900 -1624.512992 

#56: NPV ExtractVal ues , 10 

5 -x - 1200 -1924.512992 

6 -x - 1500 -2224.512992 

7 -x - 1800 -2524.512992 

8 -x - 2100 -2824.512992 

9 -x - 2400 -3124.512992 

s 10 -x - 2700 -3424.512992 . 

Nominal Interest Rate: 10% 

Net Present Value: -11319.22133 

#57: 

Equivalent Periodic Value: -1842.151146 

Net Future Value: -29359.14499 

If you have any quesitons or suggestions for improvement, please feel free to 
send emai 1 . 

Financial Mathematics IV will cover the B1 ack-Schol es or European Options Model 
and American Call Options using the bivariate normal distribution. 

[1] Financial Mathematics I (Bond, Price and Yield), DNL#52 

[2] Financial Mathematics II (TVM and Amortization), DNL#55 
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The following contribution enable simulating random games using DERIVE-functions. The 
very special parts are the graphic realisation of coins, dice, wheel of fortune and roulette. 
Studying Lorenz' files is very useful for those among you who want to learn how to produce 
graphic representations. More can be found on his very recommendable website 

http://lzkopp.gmxhome.de/index.htm 

We will present more examples in one of the future DNLs. Josef 

Four Graphic Simulations of Random Games 

Lorenz Kopp 

#2 - #14: Utility-file to produce the numbers 0 to 99 

(Prepared to save it as a mth-file. You find it saved as numbut . mth among the files) 

1. "Throw Coins" 

rd(k_): Random- number 1 ... k_ 

kd: Constant with last value of rd(k_) 

coi ngraph(c_, r_, i_): Circle marked with "1" or "0" 

Settings of 2-D-plot window for the following three examples: 

Tile vertically, reset plot range, zoom out once, reset aspect ratio, cross position (7 / 7), 
cross as center (-1 < x, y < 15), simplify before plotting, points connect, points small, grids off. 

Next 3 expressions: 

1. One "coin" center (8/8), r = 4, 

2. One coin, center (by random function) 4/7 < x,y < 100/7, r = 1/2 (repeat plotting), 

3. 20 throws of a coin randomly distributed 


Random(O) 

coi ngraph([8 , 8] ,4,36) 


12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 



1 2 3 4 5 6 7 S 9 10 11 12 


coi ngraph 


ran(4,100) ran(4,100) 


1 
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2. "Wheel of Fortune" 

rotf i g (f_ , d_) : Rotates a vector f_ of polygons, center (0/0), angle d_ 

wheel g raph ( k_ , r_) : Circle with k_ sectors, radius r_, center (0/0), a marking triangle on the top, 
area of circles shaded from 0 to 0.7r_, step 0. lr_ 

(The third element VECT0R(crcarO . . . ) only seems to imitate a rotation, it may be deleted.) 
rotnum(k_, r_): Plots the numbers 1 to k_ equidistant on a circle, number kd on the highest point 
rotwheel (k , r_): Executes wheelgraph, random-number rd and rotnum 

2D-plot-settings for the following example: 

(r_ < 1 to avoid plotting the line y = kd (1 < kd < k_) in the plot range) 

Axes: lines, labels, titles off, Grids: off, Points: connect, small, Approximate before plotting. 
Windows: tile vertically, plot range reset, zoom in twice (-0.8 < x, y < 0.8), reset aspect ra- 
tio. 

(The preset values in #6 can be used till k_= 30 (e.g. 31 ... 45: yu:=0.85, fh:=0.08).) 
random(O) 


rotwheel (15 ,0.7) 


rotwheel (10 , 1) 




3. "Play Dice" with 1 to 4 Dice 

pt(c_, r_): #32 points plotted as 4 circles (faster) or #33 points shaded (slower) 
al 1 pt (c_ , a_) : 7 points with order 4 3 

6 1 7 
2 5 

di ce(c_, a_, k_): Square with center c_, length 2a_ with k_ points 

rn(n_): Simplify produces a vector of n_ random numbers 1 ... 6, 

kn: Constant with last value of rn(n_) 

cent(a_): Centers of 4 squares (1.-4. quadrant) 

play(a_, n_): Plots n_ dice 

playsum(a_, n_):Sumof points added 

The same settings of 2-D-plot window as in nr. 2. 

Mark one of the functions in the next two lines and plot. 
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random(O) 

[pi ay (0.3,1) , pi ay (0.3, 2) , pi ay (0.3, 3) , playCO . 3 ,4) ] 

[ pi aysum (0.3,1) , pi ay sum (0 . 3 , 2) , pi ay sum (0 . 3 , 3) , pi aysum(0 . 3 , 4) ] 





• 

□ 

• 


• # 


• 








• 


w w 


# 


HP 


pi aysum(0 . 3 , 2) pi aysum(0 . 3 ,4) 


4. "Roulette" 

list: List of roulette-numbers 

rotl i st ( r_) : Numbers plotted equidistant on a circle 

noi r(r_): 37 radial lines, black sectors marked with 8 shorter lines 

roulette (r_) plots 3 circles, radial lines and numbers (roulette wheel) 

pell et( r_) plots dependent on random number a shaded circle under one number 
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Plot first roul ette (slow), then pel 1 et (fast). Delete last plot (only pel 1 et) and repeat plotting pel- 
let. 

Set yu: =0.85, f h : =0 . 08. 

[ roul ette(0 . 7) , pel 1 et(0 . 7) ] 



Derive- and CAS-TI-User Forum 2 

continued from page 5 


SHAPE (a, d := {}, n, m, 1) : = 

Prog 

If d = {} 

RETURN DIMS (a) 

If a = [] 

RETURN "REPARTITION: err: Can't repartiton EMPTY" 

If -i VECTOR? (a) 

RETURN REPARTITION ( [a] , d) 

If -i VECTOR? (d) 

RETURN REPARTITION (a, [d] ) 

If d = [] 

RETURN (RAVEL (a) ) &#8595; 1 
n := & # 8 7 1 9 ; (d) 
m := DIM (a) 

1 := DIM (d) 

a := (ITERATE (APPEND (a, x) , x, a, n/m) ) &#8595; [1, ..., n] 

(ITERATE ( [PARTITION (x, d&#8595;y), y - 1], [x, y] , [a, 1], 1-1)) 

& # 8 5 9 5 ; 1 

This now works like the "reshape" (rho) function in APL, monadically and dyadically. 


-drl 
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Lunsford 21.2.2006 

One more refinement to make it recognize strings: 

SHAPE (a, d := {}, n, m, 1) : = 

Prog 

if d = {} 

If STRING? (a) 

RETURN [DIM (a)] 

RETURN DIMS (a) 

If a = [] 

RETURN "SHAPE: err: Can't reshape EMPTY" 

If NOT VECTOR? (a) 

RETURN SHAPE ([a], d) 

If NOT VECTOR? (d) 

RETURN SHAPE (a, [d] ) 

If d = [] 

RETURN (RAVEL (a)) SUB 1 
n := PRODUCT (d) 
m := DIM (a) 

1 := DIM (d) 

a := (ITERATE (APPEND (a, x) , x, a, n/m) ) SUB [1, n] 

(ITERATE ( [PARTITION (x, d SUB y), y - 1], [x, y] , [a, 1], 1-1)) 

SUB 1 

Set Operations 

Dr. Lunsford, 17.2.06 08:39 

Derive 6.1 gives {m} A {n} as {[e A N, M]}. You can extract the member with FIRST 

FIRST ( {m} A { n } ) is [e A N, M] 

So FIRST({m} A {n}) SUB 1 is e A n etc. 

This behavior seems consistent, if you add elements you get for example 
{M,N,P,Q} A {A,B,C,D} returns 

{ [e A c, N] , [e A c, m] , [e A c, p] , [e A c, q] , [e A D, N] , [e A D, m] , [e A D, p] , 
[e A D, q] , [e A B, N] , [e A B, m] , [e A B, p] , [e A B, q] , [e A A, N] , [e A A, m] , [e A A, 
p] , [e A A, q] } 

What on earth is this about? Note the odd order of symbols, and their capitalization. This behavior is 
not documented with the product. 

If you think of A A B = e A (B log A) we can write in Derive 

[e A N, M} as e A {N, log M} 

so this in a sense is the explanation of {M} A {N}. 

I assume all this is well known? 

Any pointers to online resources would be appreciated. I notice this list is not very active. Is there a 
compiled bug list? Idiom list? What a frustrating product. 


-drl 
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Valeriu Anisiu 17.2.06, 15:37 

If A,B are sets, A A B is not defined in Derive, so it should be returned unevaluated, (only A A n makes 
sense, where n is a positive integer). The fact that A A B returns some garbage is a bug but not a very 
serious one. 

V. Anisiu 

Confusion with Complex Numbers - imaginary part zero!! 

Heinrich Ludwig 28.2.2006 

Sehr geehrter Herr Bohm, 

ich bin auf einen sehr subtilen Fehler in DERIVE 6.1 gestoften, der wegen seines gehassig versteck- 
ten Auftretens vielleicht noch nicht bekannt ist. 

Bemerkenswerterweise tritt er in DERIVE 5 nicht auf. DERIVE 6.1 lauft bei mir auf Windows 98, Built 
4.10.2222A. 

Bitte betrachten Sie folgenden DERIVE-Text. Da mth-Dateien dankenswerterweise noch ASCII-Text 
sind, habe ich der Einfachheit halber den Code hier eingefugt. Ich bitte Sie, ihn fur das Weitere in eine 
mth-Datei zu kopieren. 

F (x, m, p) : =PROG (m:=x/ (-3 + 0*#i) , p : =m A 2 , RETURN ( [m,p] ) ) 

F ( [ [-3, 6] , [-6,0]]) 

; Simp (#2 ) 

[[[1,-2], [2,0]], [[-3,-2], [2,?]]] 

Wenn Sie die Funktion F mit der angegebenen Matrix, also #2, auswerten, erhalten Sie zwei Matrizen. 
Die erste ist, wie man dem Code entnimmt, der (-3)-te Teil des Funktions-Arguments. Soweit ist alles 
in Ordnung. Quadriert man diese Matrix wiederum, so stellt sich der Fehler ein. Er besteht darin, dass 
-3+0 #i intern nicht gleich -3 zu sein scheint. 

Allerdings kann man durch keine der ublichen Maftnahmen (DISPLAY, Vergleich mit -3 ergibt true, 
usw.) den Unterschied erkennen. Auch die Matrix m schaut korrekt aus. Erst beim Quadrieren zeigt 
sich, daft m doch nicht gleich [1, -2; 2, 0] ist, obwohl es nach auften den Anschein hat und Tests der 
Koeffizienten mit INTEGER?( ) true ergeben. Auch INTEGER?(3+0*#i) wird vernunftigerweise zu true 
vereinfacht. 

Ich bin auf den Fehler gestoften, als ich die von Herrn Prof. Wiesenbauer veroffentlichte Funktion 
EXTENDED_GCD verwendet habe. Sie enthalt die Zeile 

q_ := ROUND (RE (q_) ) + ROUND ( IM (q_) ) *#i 

Dabei kommt es durchaus vor, dass der Imaginarteil von q_ Null ist. 

Verwendet man den Output der Funktion weiter, dann bekommt man das oben geschilderte Problem, 
jedenfalls mit DERIVE 6.1. Naturlich hat Prof. Wiesenbauer vollkommen korrekt programmiert und zu 
DERIVE-5-Zeiten hat sich der Code auch noch richtig verhalten. 

Heinrich Ludwig came across a strange fact, that working with a complex number having the imagi- 
nary part 0 (i.e. a + O i) can lead to serious problems. Internally a and a + 0-i seem to be different as 
one can see following the matrix-operation given above (matrix m divided by -3+0 i and then squared). 
Heinrich found this bug? first when using Johann Wiesenbauer's EXTENDED GCD which can show 
0 as imaginary part. Using the output of this function leads immediately to the problem. In DERIVE 5 
one cannot find this strange behaviour. 


p34 


Derive- and CAS-TI-User Forum 


D— N— L#61 


/\s Mr Ludwig refers to Johann 's function I sent this mail to him and as usual he immediately resent an 
answer: 

Johann Wiesenbauer 2.3.2006 

Lieber Josef! 

Leider muss ich Herrn Ludwig recht geben, es liegt hier ein offensichtlicher Bug in der aktuellen Ver- 
sion von Derive vor. In seiner einfachsten Form kann man ihn erkennen an der Funktion h(x), die fol- 
gendermalien definiert ist: 

h(x) := PROG(x := x-(l + 0-#i), [x, x-x]) 

Damit ist z.B. h([l, 2]) = [ [ l , 2], 5] 

was durchaus korrekt ist, aber eben h([l, 0]) = [ [ l , 0], ?] 

d.h. das Problem ist, dass hier zwei Nullen aufeinanderstossen, die intern offenbar als verschieden 
behandelt werden. 

Ich kann daher nur empfehlen, diesen Bug in der Derive 6.1 umgehend an Albert Rich weiterzuleiten, 
damit er im nachsten Patch behoben wird. 

Johann Wiesenbauer 2.3.2006, 5 hours later 

Hier noch ein einfacheres Beispiel fur den Bug, welches zeigt, dass das Ganze mit Vektoren und Mat- 
rizen gar nichts zu tun hat. 

g(x) := PROG(x :* 1 + 0#i, [x, x-x]) 


Es ist 

g (5) 

= [5, 25] 

Aber 

g (0) 

= [0, ?] 


Bis zum nachsten Patch, in dem das hoffentlich behoben sein wird, fallt mir auch keine bessere vor- 
laufige "Reparatur" als 

g(x) := PROG(x :* 1 + 0#i, x:=if (x=0, 0, x) , [x, x-x]) 

ein. 

Johann sent two mails within 5 hours confirming Heinrich Ludwig's findings. This is a bug and it 
should be fixed in the next update of DERIVE. The problem is that two kinds of zeros are appearing. 
The bug has nothing to do with vectors and matrices. This is the content of Johann's 2 nd mail giving a 
possibility to “repair” this deficiency for his example of a multiplication by 1+0 -i (see above). 

Heinrich Ludwig 3.3.2006, 12:15 

Sehr geehrter Herr Bohm, 

vielen Dank, dass Sie so schnell auf meine jungste Fehlermeldung reagiert haben. Das ermutigt mich, 
Ihnen von einem weiteren Bugs in DERIVE zu berichten. Ich werde Sie noch liber einen weiteren 
Fehler unterrichten, brauche aber noch etwas Zeit zur genauen Analyse. 

Der neue Fehler hangt mit dem vor kurzem geschilderten zusammen, ist aber nicht mit ihm identisch, 
da er sowohl in DERIVE 5.06 als auch in DERIVE 6.1 auftritt. Er hat offenbar damit zu tun, dass 
DERIVE Schwierigkeiten hat, komplexe Zahlen, deren Imaginarteil gleich Null ist mit reellen Zahlen zu 
identifizieren. Ich mochte nicht sagen, dass DERIVE das nicht konne Oder tue, sondern nur, dass 
diese Identifikation verzogert geschieht, was beim Programmieren zu Schwierigkeiten fiihrt, die zu 
umgehen unnotige "Workarounds" erfordert. 
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DERIVE erkennt nicht von Haus aus, dass z*CONJ(z) stets reell ist, sondern fuhrt die Konjugation 
und anschlieliend die Multiplikation aus. Dabei ergibt sich (aus Sicht von DERIVE zufallig), dass der 
Imaginarteil Null ist. Ubergibt man das Ergebnis einer Funktion, die prinzipiell nur reellwertige Argu- 
mente akzeptiert, wie MAX, MIN, GCD, LCM etc., so sollte dies fur DERIVE Anlass sein, das uberge- 
bene Argument versuchsweise zu einer reellen Zahl zu vereinfachen. Das tut DERIVE merkwurdiger- 
weise 

nicht konsequent. Sie finden das dargestellt in einer mth-Datei im Anhang dieses E-Mails. Interessant 
ist, dass die von den Funktionen a und b erzeugten Matrizen bei weiterer Vereinfachung (CTRL-B) 
richtig sind, d.h. dass nicht zu Ende vereinfacht wurde. Hierin unterscheidet sich der heute gemeldete 
Bug von dem vorgestern gemeldeten wesentlich. 

Naturlich lasst sich das Fehlverhalten leicht vermeiden, wenn man das Betragsquadrat der komplexen 
Zahl mit RE(-) kapselt Oder gleich stattdessen RE(z) A 2 + IM(z) A 2 schreibt. Dennoch bleibt ein Fehler 
ein Fehler. 

Ich mochte auch gerne "Positives" beitragen. Mir gefallen Ihre Tipps und Tricks recht gut, die Sie ver- 
offentlichen, und hatte da ein bisschen was, was vielleicht dazupassen wurde. Das mochte ich Ihnen 
aber in getrennter Post zuschicken, sonst wird diese Nachricht zu lang. 

Mit freundlichem Gruli 
Heinrich Ludwig 

FEHLER. MTH: Fehl erdemonstrati on . 
f(z) :z MAX(1, z-CONKz)) 
g(z) := GCD(1 , z-CONKz)) 
h(z) := GCD(0 , z-CONKz)) 

a := VECTOR (VECTOR(f(x + L-y), x, -2, 2), y, -2, 2) 

b := VECTOR (VECTOR(g(x + L-y), x, -2, 2), y, -2, 2) 

c := VECTOR (VECTOR(h(x + L-y), x, -2, 2), y, -2, 2) 

a 

8 5 MAX(1, 4) 5 8 

5 2 MAX(1, 1) 2 5 

MAX(1 , 4) MAX(1 , 1) MAX(1, 0) MAX(1, 1) MAX(1, 4) 

5 2 MAX(1, 1) 2 5 

8 5 MAX(1, 4) 5 8 

b 

1 1 GCD(1, 4) 1 1 

1 1 GCD(1, 1) 1 1 

GCD(1 , 4) GCD(1 , 1) GCD(1, 0) GCD(1, 1) GCD(1, 4) 

1 1 GCD(1, 1) 1 1 

1 1 GCD(1, 4) 1 1 

c 

'8 5 4 5 8' 

5 2 12 5 

4 10 14 

5 2 12 5 


8 5 4 5 8 
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In his second mail Heinrich Ludwig reports about another bug? which occurs not only in DERIVE 6.1 
but also in DERIVE 5.06. DERIVE seems to have difficulties to identify complex numbers with imagi- 
nary part 0 with real numbers. It is not that DERIVE cannot do this, but it needs an additional simplify- 
ing step which can lead to problems within a program and requires some unnecessary “workarounds”. 
So DERIVE does obviously not recognize that z * CONJ(z) always gives a real result, but first finds 
the conjugate and then performs the multiplication. This leads - by chance from DERIVE' s point of 
view - to imaginary part 0. Using the product in a function which accepts only real arguments like 
MAX; MIN; GCD, LCM etc, DERIVE should try to simplify it to a real number, but DERIVE does not. 
You can find a documentation of this behaviour in the attached MTH-file (see above). It is very inter- 
esting that the matrices produced by a and b will become correct after an additional simplify 
(CTRL+B). 

Of course one can avoid this incorrect behaviour by encapsuling the squared absolute value with RE( 
) or writes RE(z) r '2+iM(z) r '2. But in my opinion a bug remains a bug. 

Regards 

Heinrich Ludwig 3.3.2006, 14:00 

Sehr geehrter Herr Bohm, 

nachdem ich Ihnen in kurzer Folge geschrieben habe, entschuldigen Sie bitte, dass ich's nochmal 
tue. Ich habe als Lehrer in Bayern gerade Ferien und das ist die einzige Zeit, in der man endlich ein- 
mal Zeit hat, Dinge nicht nur oberflachlich abzutun. Ich mochte Sie keineswegs in irgendeiner Weise 
drangen, weil ich so kurz aufeinanderfolgend schreibe. 

Wussten Sie schon, dass 5374978561 - 2403763488-V5 = 0 ? DERIVE 6.1 ist gelegentlich der Mei- 
nung, auch im Vereinfachungsmodus "Exact". Das Verhalten trat in der Version 5 noch nicht auf, wie 
ich unten darlege, deshalb mochte ich es als Fehlverhalten einstufen. 

Wenn Sie in der anhangenden Datei Fehler3.mth nachsehen, erkennen Sie in der Funktion GGT den 
Euklidischen Algorithmus. Zu Analysezwecken ist die Display-Anweisung eingefugt. Wendet man den 
Algorithmus auf zwei nicht kommensurable Zahlen an, z.B. 1 und sqrt(5), dann sollte die Abbruchbe- 
dingung b=0 nie erfiillt sein. Mit DERIVE 6.1 terminiert GGT aber nach endlich vielen Schritten, weil 
aus unbekannten Grunden die o.g. 

merkwurdige Gleichheit eintritt. Die Erscheinung ist nicht singular, man kann leicht auch andere Wur- 
zeln nehmen. Sicherlich ist die o.g. Differenz eine recht gute Approximation an Null (9.3E-11), aber 
DERIVE sollte doch "exact" rechnen! Im Modus "mixed" ist an DERIVES Verhalten nichts auszuset- 
zen. Interessant ist: Wenn man die o.g. Differenz alleine, also ohne den Kontext der Funktion GGT, 
eingibt und dann vereinfacht (Strg-B), wird sie nicht zu null, genau wie es sein soil. 

Da ich schon einmal die Erfahrung gemacht habe, dass sich von Version 5 auf 6 allerhand Internes 
verandert hat, habe ich auch Version 5 das Problem vorgelegt. Zunachst erkennt man, dass die alte 
Version im Gegensatz zur neuen ROUND(sqrt(5)) im Modus "exact" nicht weiter vereinfacht. Insofern 
ist die neue Version eine Verbesserung. Urn Version 5 die Funktion GGT sinngemalJ ausfiihren zu 
lassen, muss man zum Runden ROUND(APPROX(a/b, n) setzen. (Dies ist in GGT2 ausgefiihrt.) Nun 
zeigt sich 

Erstaunliches: Beide Versionen verhalten sich 17 Schleifendurchlaufe lang gleich. Wahrend Version 
6.1 nach der 17. abbricht, rechnet Version 5 weiter und kommt nur zum Absturz, weil irgendwann a/b 
bei der eingestellten Stellenzahl zu "+-"inf approximiert wird. Also enthalt Version 5 den genannten 
Fehler nicht. 

Darf ich ehrlich sein? Ich neige jetzt immer mehr dazu, mich wieder an Version 5 zu halten. Sollte ich 
einmal das nutzliche Degug-Feature DISPLAY brauchen, so werde ich fur zwischendrin ab und zu die 
Sechs-Einser nehmen. 


Schade. 
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Heinrich Ludwig's 3 rd report: In his mail he points out that DERIVE sometimes has the opinion that 
5374978561 - 2403763488 ^5 = 0 even in Exact Mode. In the attached file FEHLER3.MTH you will 
recognize the Euklid algorithm in function GGT. For analysing purpose I inserted the DISPLA Y instruc- 
tion. Applying the GGT on two incommensurable numbers as 1 and sqrt(5) then the process should 
never end. DERIVE 6. 1 terminates GGT after some steps. DERIVE 5 produces an infinite loop which 
is correct. 


Surely 9.3E-1 1 is a fine approximation for zero but I expect DERIVE to work “EXACT” in Exact Mode. 
In Mixed Mode everything is ok. Interesting detail simplifying the above difference (CTRL+B) does not 
result in 0. 

Version 5 does not simplify ROUND(sqrt(5)) in Exact Mode - this is an improvement in version 6. For 
performing GGT in version 5 one has to set ROUND(APPROX(a/b,n)) - see in GGT2. We can ob- 
serve the following: 

Both versions show the same behaviour for the first 17 loops. Version 6.1 stops then but version 5 
goes on calculating and crashes at a certain moment because a/b will be simplified to “+-inf”. So De- 
rive 5 works exact. 


Precision := Exact 
Preci si onDi gi ts := 10 
Notation := Decimal 
GGT(a, b, q, r) := 

Loop 

If b = 0 
RETURN a 
q := R0UND(a/b) 
r := a - q • b 
a := b 
b := r 

DISPLAY([b, APPR0X(b, 27)]) 

GGT (1 , V5) 

[ 31622993-^5 - 70711162 , 0 . 000000007071019424 ] 

[299537289 - 133957148 -^ 5 , 0 . 000000001669241254 ] 

[ 567451585 -JS - 1268860318 , 0 . 0000000003940544068 ] 

[5374978561 - 2403763488 -^ 5 , 0 . 00000000009302362684 ] 

#7: GGT (1, VS) = 567451585-^5 - 1268860318 

#8: GGT 2(1, ^5) = [5374978561 - 2403763488-^5, 567451585-^5 - 1268860318] 


GGT2(a, b, q, r, i := 0) := 

Loop 
i :+ 1 
If b = 0 
RETURN a 

q := ROUND (APPROX (a/b, 27)) 
r := a - q • b 
a := b 
b := r 
If i = 17 

RETURN [b, a] 


Heinrich Ludwig 27.3.2006 

Ich habe wieder einen Derive-Bug zu melden. Er steht moglicherweise in Zusammenhang mit denen, 
uber die ich jiingst berichtet habe, aber es muss nicht so sein. Leider bin ich jetzt zu spat dran, dass 
im DNL#61 noch daruber berichten werden konnte. Ich bin erst gestern von meinem Kollegen Wolf- 
gang Jager - einem langjahrigen treuen Anhanger von Derive ubrigens, auf ihn geht auch die Ver- 
wendung von Derive an unserer Schule zuriick - darauf hingewiesen worden. 
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COS(ACOS(x)) wird zu x vereinfacht, wenn man den Definitionsbereich nicht einschrankt. Das ist in 
komplexer Rechnung richtig und Derive macht gute Arbeit. Wenn man aber mit Derive 5 den Term 
COS(ACOS(x)) in einem 2D-Fenster sich graphisch darstellen lasst, dann wird derTerm nurfiirx aus 
[-1; +oo[ gezeichnet. Derive 6 zeichnet dagegen richtig. 

Nun wird es kein Patch mehr fur Derive 5 geben, nachdem man fur Derive 6 schon so lange auf Pat- 
ches wartet. Aber meine Fehlermeldung mag die Programmierer von Derive darauf hinweisen, noch 
einmal genau nachzuschauen, ob in Derive 6 der Fehler wirklich behoben ist, Oder ob nur ein Sonder- 
fall vorliegt, der diesmal das Fehlverhalten kaschiert. Wenn Sie sich moglicherweise erinnern, so ha- 
ben die Bugs, die ich gemeldet habe, alle damit zu tun, dass a + 0*#i von Derive nicht zuverlassig als 
reelle Zahl erkannt und weiterverarbeitet wird. Genau darum scheint es sich auch hier zu handeln: 
Schreibt man RE(COS(ACOS(x))), so tritt der Fehler nicht mehr auf und Derive 5 zeichnet korrekt. 

Heinrich Ludwig points out that there is another problem which might be connected with the bugs 
mentioned earlier: 

COS(ACOS(x)) is simplified to x if the range is not restricted. This is ok for calculations with complex 
numbers and DERIVE does a good job. Plotting the graph with DERIVE 5 the graph only appears for 
[-1,+ oc[ : the DERIVE 6 plot is correct. The developers of DERIVE should investigate (for a hopefully 
coming Update) if the bug has really been fixed or if there is only a special case which hides the bug. 
I believe that this behaviour again is based on the fact that DERIVE does not recognize in many cases 
that a + 0*#i is a real number. Writing RE(COS(ACOS(x)) solves the problem in DERIVE 5 and leads 
to a correct plot, too. 

Manuel Garcia Mateos 22.3.2006 

Hello Derivians, 

I tried to solve the following complex equation of third degree 2z A 3-(3+15#i)z A 2+(-45+9#i)z+54#i=0 
One solution should be z=3#i (the others are 3+3#i and 3/2+3/2#i). 

I don't get these solutions. 

After a substition of z by 3#i (subst(eq,z,3#i)) Derive says "false". 

Where is the error? 

Greetings, 

Manuel 

Aleksey Tetyorko 

Hi! 

I do not understand your troubles. My Derive 4.1 1 solves: 

2*z A 3-(3+15*#i )*z A 2+(-45+9*#i )*z+54*#i =0 

; Sol ve(#l) 
z=3*#i 

; Sol ve(#l) 
z=3+3*#i 

; Sol ve(#l) 
z=-3/2+3*#i/2 


Aleksey 
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Johann Wiesenbauer 23.3.2006 

Hi Manuel, 

It is true, I had also troubles when trying to solve your equation algebraically in exact mode. To be 
more precise, I got a screen full of nested square roots, but Derive doesn't see how to simplify them. 
It can see though that 3#i is a solution if you check it like this 

3 2 

w := 2-z - (3 + 15-O-z + (-45 + 9 -l)-z + 54 - ZL 

SUBST(w, z, 3 ■ t.) = 0 

With "a little help" it can also compute the exact solutions using the following small program 

mysolve(w) := [1, l] •S0LUTI0NS([RE(SUBST(w, z, u + v-t.)), 

IM(SUBST(w, z, u + v-l))], [u, v] , Real)' 

3 3 • L ' 

mysolve(w) = 3-L, 3 + 3 • L , - + 

2 2 . 

Hope this helps. 

Cheers, Johann 

Stefan Welke 23.3.2006 

Hello Manuel and Johann, 

in contradiction to Johann's observation my DERIVE 6 solves this equation: 

Just try 

3 2 

S0LVE(2-z - (3 + 15-O-z + (-45 + 9 -l)-z + 54-l, z) 

The result is very clumsy and contains many roots etc. Then apply 
Simplify->E3asic from the menu bar and you will immediately get: 

3 3 • L 

z = - + vz = 3 + 3- lvz = 3- l 

2 2 

Greetings, 

Stefan Welke 

Johann Wiesenbauer 23.3.2006 

Hi Stefan, 

While it is true that Derive 6.1 can surprisingly simplify its first result when manually forced to do so, 
I'm not sure if this a "contradiction" to what I said. At any rate, this "2-step simplification" is a little bit 
unsatisfactory. In a way Derive behaves like a pupil, who is unable to produce the expected result in 
the first try, but can do so after being given the hint that a further simplification is possible. 

Cheers, Johann 


Wim de Jong 24.3.2006 

Hi all, 

It is indeed surprising that the complicated algebraic solution is not directly simplified. However if you 
Simplify (Basic) it you get the required answer. More directly, Approximation of the expression 
Solve(2z A 3-(3+15i)z A 2+(-45+9i)z+54i=0,z) does the trick. 

As to the substitution of 3i for z, Simplify>Subexpression Substitution yields 0=0. 

Cheers, 

Wim de Jong 
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Two strange limits 

Jim FitzSimons 

1 + COS(x) 

This limit is wrong: i™ =1 

X-KO 1 - COS (x) 

David Sjostrand 

Derive doesn't seem to be able to compute the limit LIM(f(x + 1)/f(x), x, inf), with f(x) := 2 A x + 3 A x 


X X 

f(x) := 2 + 3 

(x + l)/x 10 (x + l)/x 10 

f(x +1) 1 -LN(3) + 2 • LN(2) 

1 ill = lin 

x-k» f (x) X-.0+ 1/x 10 1/x 10 

3 • LN(3) + 2 • LN(2) 

1/x 10 1/x 10 

f(x + 1) 2-3 • LN(3) +3-2 -LN(2) 

lira = lira 

x->— <o f(x) x->0+ 1/x 10 1/x 10 

3 • LN(3) + 2 -LN(2) 


Let's look how the Tl deals with the problems from above: 



Jan Vermeylen 1 1 .3.2006 

Hello Josef, 

I have a question of one of my science colleagues. 

He has a set of some thousand quadruples (t,x,y,z) which show the position (x,y,z) of a particle at a 
certain time t. 

Is it possible in Derive to produce a 3D graph of the path of this particle with an animation of the mov- 
ing particle? Or should I use DPGraph? 

Thank you for giving me advice. 

Best regards, 

Jan 
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DNL: Good question, Jan. I tried to introduce a slider bar for the subscript of a point matrix to 
address one point after the other for being plotted but I failed. Obviously the slider bar cannot 
be applied on discrete variables? 

It would be great if any of our graphic sorceres had any advice. Josef 

Last News: 

Peter Schofield sent a mail containing very valueable advice - see next DNL 


An InterConnectivity-Communication problem between DERIVE and the Tl 

Rainer Gerdes from Germany asked for advice because he failed trying to transfer a matrix 
equation from the Tl to DERIVE. This might be of interest for some of you. See a summary of 
his report, and how to overcome this deficiency: 

I enter the equation [l,2;4,5]*[x;y] = [l,2] on the Voyage200, save it as a textfile, connect the 
device with my PC, launch DERIVE and import the file. After some seconds the error message “Syn- 
tax Error” appears. What am I doing wrong? 

Regards 

Rainer Gerdes 


DNL: Dear Mr. Gerdes, 

You are doing right, InterConnectivita does wrong - that's the bad news, but your matrix calculation 
can be transferred - but with some tricks, that's the good news. There is something to be improved 
by the software smithies. 

One way is to edit the matrices in the text editor of the Tl as lists. Your equation should then read as 
follows (first attempt): 

{{1,2}, {4,5}}*{{x}, { y } } = { l , 2 } (which delivers the result false because of the non suiting 
dimensions, the next line will do better:) 

{{1,2}, {4,5}}*{{x}, {y}} = {{l}, {2}} will result in simultaneous equations in DERIVE- 
environment. 

Interestingly enough one can leave the matrices (and the product but not the whole equation) with 
brackets and it will be transferred correctly. So you can also transfer without any problem the following 
forms of the equation: 


[1 , 2 ; 4 , 5 ] * [x;y]={l,2} and [ 1 , 2 ; 4 , 5 ] * [x; y] ={ { 1 } , { 2 } } . 

And there is another very "perfidious" trick: 

In the text editor remove the leading "C"s (which declare the lines as executable commands => leading 
to commands in DERIVE) converting the lines to simple comments. This will be very welcome by 
DERIVE appearing as a comment in a textbox. In DERIVE copy and paste this "fictitious comment" 
into the edit line, ENTER, and .... it is done. 

Best regards and much luck 
Josef 
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Susan Jackman from England asked how to produce the DERIVE-Snail (Cover of the manual). As you 
might remember we had this snail in the last DNL. I sent the file to Susan, but she had problems to 
bring the snail alive with DERIVE 5 (no slider bars). Finally it lucked but in the meanwhile one of her 
students had made his/her own tries. Susan wrote: 

Josef, 

Thank you. I now have opened your snail and I'm very impressed. I am attaching my seashell which is 
based on a much simpler function. The idea came from a student who found something like it (done 
on Mathematica) on www.wolfram.com However I spent some time changing values of the parame- 
ters and changing the scale. 

Best wishes, 

Susan 


v 't 
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If anybody of you has other animals for our DERIVE - zoo, please send your pictures, Josef 

On the Titbits: 

We have no Titbits in this issue, because the journal is full and Johann has a lot to do 
on the University. He promised the next Titbits for DNL#62. Come and meet Johann 
in Dresden. He will enjoy any talk about DERIVE and his column. Josef 
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Submitted and Accepted Lectures and Workshops for the 7 th International 
DERIVE and TI-CAS-Conference in the Frame of DES -TIME 2006 Dresden 

(in order of submitting the papers) 


The Lectures: 

Karsten Schmidt, D, 

Alex Lobregt, NL, 

Dennis Pence, USA 

Michel Beaudin, CAN 

Jose Luis Galan, SPA 
Jose Luis Galan, SPA 
Josef Bohm, A 
Giora Mann, ISR 
Lisa Townsley, USA 
Gabriel Aguilera, SPA 
Ludwig Paditz, D 
Gilles Picard, CAN 
David Jeffrey, CAN 
Steven Schonefeld, USA 
Jean-Jacques Dahan, 
Max Schrofel, D 
Gunter Alfanz, A 
Eugenio Roanes, SPA 
Renee Gossez, BEL 
Gilles Picard, CAN 

Stephan Joubert, RSA 

Josef Lechner, A 

David Sjostrand,SWE 
Roman Hasek, CZ 
Walter Wegcheider, A 


Working with Random Variables and Random Numbers in DERIVE 
Signal Processing using a CAS 

Activities with Implicit Functions and Implicit Differentiation 
the TI-89A/200 

Theory versus Technology: don't look for a competition, 
look for collaboration 

LINEJNTEGRALS.MTH: Solving Problems of Line Integrals using Derive 
MULTIPLEJNTEGRALS.MTH: Solving Problems of Multiple Integrals 
Background Pictures as a Stimulating Means for Math Teaching 
Changing the Viewing Angle on a Conic Section 
DERIVE, a Mathematical Assistant 

GRAPH_ALGORITHMS.MTH: Graph algorithms using display step in D6 

The Rank of a Matrix with Parameters and the Solution of a Linear System 

Using the Voyage 200 in the Classroom: Surprising Results 

Flying with DERIVE 

Amorous Bugs and Pursuit Problems 

FR Vizualizing Solutions of Des and Systems of 2nd order with Cabri2+ 
Cubic Spline Functions, Matrices and DERIVE 

Der Einsatz von CBL/CBR im facherubergreifenden, naturw. Unterricht 

Eight Wishes about Computer Algebra Systems 

Using DERIVE 6 to find equation and to visualize a locus in 3D space 

Numerical Methods with the V200: To Function or to Program, 
that's the quest 

The Mathematics and Industrial Applications associated with the 
singing wine glass using a CAS 

See it complex - make it easy, Experimental working using 
complex numbers 

Interactive Investigations with Slider Bars in DERIVE 6 
Geometry with DERIVE 

Solving Problems of Spherical Trigonometry with the Help of 
Computer Algebra 
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Johann Wiesenbauer, A 
Pavel Pech, CZ 
Tania Koller, A 
Temple Fay, RSA 
Wolfgang Schlichthorn, D 
Quay van der Hoff, RSA 
Anna Salvadori, I 
Leon Magiera, POL 
Peter van der Velden, NL 


Elliptic Curve Cryptography with DERIVE 
Discovering and Proving Geometric Inequalities by CAS 
Motivierender Mathematikunterricht mit CAS - ein Erfahrungsbericht 
A Quadratic Spring Model 

Das Rechnen mit GrolJen in Anwendungsaufgaben 
Real-Life Examples of Limit Cycles 
Mathematics & Reality 

Interactive Work with CAS when Solving Problems on Physics 
Didactical Principles of Integrated Learning Math with CAS 


The Workshops: 

Rene Hugelshofer, CH 
Jose Luis Galan, SPA 

Josef Bohm, A 

Lisa Townsley, USA 
Jean-Jacques Dahan, FR 

Guido Herweyers, BEL 
Matija Lokar, SLO 
Gunter Alfanz, A E 
Carl Leinbach.USA 
Eoghan MacAogain, IRL 

Johann Wiesenbauer, A 
Wade Ellis, USA 


Dynamic Algebra 

DERIVE as a Pedagogical CAS: Programming using the 
DISPLAY-function 

Slider Bars narrow the Gap between Computeralgebra and 
Dynamic Geometry 

Calculus Lab Transformation from DERIVE to DERIVE and Biology 

Vizualizing Solutions of DEs and Systems of 1st and 2nd order 
with Cabri2+ 

Applied Linear Algebra 

Selected Themes from Teaching Calculus with CAS 

Experimente fur den Einsatz von CBL/CBR im fu, naturw. Unterricht 

Using DERIVE to Teach Bioinformatics Algorithms 

Student Projects in using DERIVE to help teach topics in 
2nd level maths 

Elliptic Curve Cryptography with DERIVE 
Integrating Graphics, Tabular, and Symbolic Calculation 


There will be another approximately exciting 60 lectures and workshop presented at the 

9 th ACDCA Summer Academy. 

You are kindly invited to visit the Conference website for browsing the scientific program, 
the venue and the social program 


www.des-time- 2006 .de 






